Quantum phase transitions of the asymmetric three-leg spin tube 
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We investigate quantum phase transitions of the S = \ three-leg antiferromagnetic spin tube with 
asymmetric inter-chain (rung) exchange interactions. On the basis of the electron tube system, we 
propose a useful effective theory to give the global phase diagram of the asymmetric spin tube. In 
addition, using other effective theories we raise the reliability of the phase diagram. The density- 
matrix renormalization-group and the numerical diagonalization analyses show that the finite spin 
gap appears in a narrow region around the rung-symmetric line, in contrast to a recent paper by 
Nishimoto and Arikawa [Phys. Rev. B 78, 054421 (2008)]. The numerical calculations indicate 
that this global phase diagram obtained by use of the effective theories is qualitatively correct. In 
the gapless phase on the phase diagram, the numerical data are fitted by a finite-size scaling in the 
c = 1 conformal field theory. We argue that all the phase transitions between the gapful and gapless 
phases belong to the Berezinskii-Kosterlitz-Thouless universality class. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Kz, 75.40.Cx 



I. INTRODUCTION 

The spin tube^^^ 10 i u i 12 i 13 i 14 i 15 i 16 i.e., the 
spin ladder with the periodic boundary condition along 
the inter-chain (rung) direction, is one of interesting mag- 
netic systems which is expected to exhibit some exotic 
phenomena due to special topology such as the carbon 
nanotube and the magnetic frustration. Actually, several 
theoretical works 6 ^ 9 ' 11 ' 12 have shown that the S = | 
three-leg antiferromagnetic spin tube has a spin gap, in 
contrast to the corresponding three-leg spin ladder with 
the open boundary condition along the rung. Moreover, 
it has been predicted in the field theoretical metho d 13 : 14 
that both the gapless and gapful vector-chiral long-range 
orders emerge in certain parameter regions of the three- 
leg tube in a magnetic field. Recently some candidates 
for the spin nanotubes have been synthesized; a three- 
leg tube [(CuCl2tachH) 3 Cl]Cl 2 (Ref.Q]), a nine-leg tube 
Na 2 V 3 7 (Ref . 0) , and a four-leg tube Cu 2 Cl 4 • D 8 C 4 S0 2 
(Refs. H and0). It is therefore expected that novel, in- 
triguing phenomena will be detected in these materials. 

Let us here focus on the S — \ three-leg antiferro- 
magnetic spin tube that is the simplest tube with geo- 
metrical frustration. Since the unit cell consists of three 
spins with S = h in this tube, the Lieb-Schultz-Mattis 
theorem^ suggests that the spin gap must be accompa- 
nied with at least doubly-degenerate ground states. In 
fact, previous numerical analyse s 6 : 11 ' 12 have confirmed 
such doubly-degenerate 5 = ground states due to 
the spontaneous breaking of the translational symme- 
try along the leg direction. The ground states have a 
valence-bond type (superposition of spin-singlet pairs) 
order j 6 -^ Here if one of the three rung coupling constants 
is changed, the following two models are reproduced as 



limiting cases: the three-leg spin ladder and the decou- 
pled system of a single chain and a two- leg ladder. These 
two systems are believed to possess a gapless excitation. 
Therefore, the S = \ three-leg system, where one of the 
three rung couplings is varied, would undergo a quan- 
tum phase transition from the gapless state to the gapful 
symmetry-broken one. A recent numerical workii has 
suggested that the gapful phase is extended to a finite 
(although narrow) region when the rung-coupling asym- 
metry is introduced. Unfortunately, however, the feature 
of the transition was not so clarified because the system 
size used in Ref. [ll| was too small. On the other hand, 
a recent density-matrix renormalization-group (DMRG) 
approach^ assuming a special power-law form of the 
finite-size correction, has concluded that the transition is 
of the first-order and the system is always gapless except 
for the symmetric three-leg spin tube. Such a discontin- 
uous transition, however, has not been reported so far in 
any realistic systems. In addition, from the viewpoint of 
the effective theory, the occurrence of such a transition 
must require a highly fine tuning of parameters. Thus the 
phase diagram and the critical properties of the quantum 
phase transitions in the three-leg spin systems are still 
controversial. 

Motivated by the above situation, in this paper, we 
study the wide ground-state phase diagram and the uni- 
versality classes of the quantum phase transitions in 
S = i three-leg antiferromagnetic spin tube with the 
rung-coupling asymmetry. We first propose a simple ef- 
fective theory to explain the quantum phase transitions 
between gapful and gapless phases on the basis of the 
Hubbard model on the tube lattice and the non-Abelian 
bosonization . 18 ' 19 : 20 This effective theory enables us to 
draw a global phase diagram by counting the number of 
Fermi points. We find one gapful phase and three gap- 
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less phases in the phase diagram, and predict that the 
gapless phases are all described by a level- 1 SU(2) Wess- 
Zumino-Witten (WZW) field theory, 1 ^ 1 ^ wmcn is a 
c = 1 conformal field theory (CFT). Besides this effective 
theory, using other analytical strategies, we consider the 
strong-rung-coupling regime and the weak-rung-coupling 
regime with a strong rung distortion (i.e., asymmetry) in 
the phase diagram. In the former regime, two of three 
gapless phases are predicted by using a known effective 
theory J* 2 - In the latter regime, we prove that a finite gap- 
less phase definitely exists. 

We subsequently perform the numerical diagonaliza- 
tion and the DMRG calculation combined with some 
finite-size scaling analyses on the basis of the above ef- 
fective theories. Applying the CFT approac h 21 ! 22 ' 23 to 
our numerical data, we argue that the quantum phase 
transitions belong to the Berezinskii-Kosterlitz-Thouless 
(BKT) universality class . 24 ' 25 A numerically quantitative 
phase diagram is presented, and is consistent with that 
of the effective theories. 

This paper is organized as follows. In Sec. [Til we define 
the Hamiltonian of an asymmetric three-leg spin tube 
model. In Sec. IIII| we draw a qualitative but global 
ground-state phase diagram by means of an effective the- 
ory based on the half-filled Hubbard model on the tube 
lattice. Employing another approach based on the non- 
Abelian bosonization, we precisely show the existence of 
a gapless phase in the weak-rung-coupling regime with 
a strong asymmetry. Moreover, we discuss the gapful 
phase in the strong-rung-coupling regime. Section IIVI is 
devoted to the numerical analyses for the spin tube. We 
plot a scaled gap calculated by the DMRG method to 
confirm the phase diagram obtained by the effective the- 
ories. We further analyze the numerical data obtained 
in the exact diagonalization on the basis of the finite-size 
scaling in the c = 1 CFT. We provide the summary and 
short discussions in Sec. [V] 



II. MODEL 

We consider the S = \ asymmetric three-leg spin tube, 
shown in Fig. [TJ described by the Hamiltonian 

3 L 

H= Ji EE^'^+i 

i=l j=l 

2 L L 

i=l j=l j=l 

where Sij is the spin-i operator and L is the length of 
the tube in the leg direction. The exchange coupling con- 
stant Ji is for the neighboring spin pairs along the legs, 
while J r and J' r are the rung coupling constants. All the 
exchange interactions are supposed to be antiferromag- 
netic (namely, positive). The ratio a — J' T I 'J r stands 
for the degree of the asymmetry of the rung couplings. 



We will vary a and J\ to investigate the quantum phase 
transitions. Throughout this paper, we fix J r to one. 




FIG. 1: Structure of the three-leg asymmetric spin tube {TJ). 



The present model includes three typical models as 
limiting cases; (a) a = 0: the three-leg spin ladder, (b) 
a = 1: the symmetric spin tube, and (c) a — > oo: the 
single chain plus rung dimers. Since the system is gapless 
in the cases (a) and (c), while gapful in the case (b), at 
least two quantum phase transitions should occur with 
increasing a from to infinity. As we already mentioned, 
the one-site translational symmetry along the leg (Sij — > 
Sij+i) is spontaneously broken in the symmetric spin 
tube at least in the strong-rung-coupling regime.— ^ 

III. EFFECTIVE THEORIES 

In this section, we study the spin tube ([TJ) by con- 
structing its low-energy effective theories. In Sec. IIII Al 
we draw a phase diagram in the whole coupling-constant 
space (a, Ji). To this end, we develop a simple effective 
theory for the spin tube from the corresponding Hub- 
bard model on the tube lattice, where the SU(2) spin- 
rotational symmetry is preserved automatically. This ef- 
fective theory allows us to find three gapless phases and 
one extended gapful phase. Next, utilizing other theoret- 
ical schemes, we carefully investigate two special regimes 
J r <C J[ <C Ji and J r » Ji , respectively, in Sees. IIII Bl 

and nro 

A. Global phase diagram derived from the 
Hubbard model on the tube lattice 

Here, we provide a systematic method to draw global 
phase diagrams of one-dimensional antiferromagnetic 
quantum spin systems. It is well-known that any S = j 
Heisenberg model is obtained from the corresponding 
half-filled Hubbard model in the limit of strong on-site 
Coulomb interactions. In one dimension, the spin config- 
urations of the low-energy states in the Heisenberg model 
qualitatively agree with those in the half-filled Hubbard 
model even with a weak Coulomb interaction. Fur- 
thermore, the phases in the weak-Coulomb regime often 
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smoothly connect to those in the strong-Coulomb regime 
in one-dimensional electron systems. Relying on these 
arguments, we construct the low-energy effective theory 
for the spin tube (Q]) from the corresponding Hubbard 
model. To discuss the wider parameter space, first we di- 
agonalize the kinetic parts of the Hubbard Hamiltonian 
including both the leg and rung hopping terms i 26 ' 27 i 28 
Then, we take account of the on-site Coulomb interac- 
tion as the perturbation, with the help of the non-Abelian 
bosonization and CFT. As one will see later, the number 
of Fermi points is essential to determine whether or not 
the Coulomb interaction open a spin gap. 

Since we consider the electron tube instead of the orig- 
inal spin tube, our results in this subsection should be 
regarded as a qualitative argument. However, the ap- 
proach from the electron tube can be applicable to the 
wide parameter space (a, Ji), in contrast to the other 
conventional methods. For example, the nonlinear sigma 
model is not usually reliable for frustrated magnets in- 
cluding the present spin tube. It is dangerous to apply 
the Abelian bosonization to any SU(2)-symmetric mag- 
nets^ Moreover, a standard non-Abelian bosonization 
method, taking into account rung couplings perturba- 
tively (we will use it in Sec. IIII B|) . is of course valid only 
in the weak-rung-coupling regime. 

Now, let us begin with the definition of the Hubbard 
model on the three-leg tube lattice. The Hamiltonian 

H = H nop + H int , (2) 

consists of the hopping part 

L 3 

Hhop — E] E] E] (^ C n+l,i,CT C «>»)<r + s *+M C n,i+l,o- C ni*><' 
n=l i=l <j=T4 

+h.c), (3) 
and the on-site interaction part 

L 3 

Hhit = U ^ ^ ra n,i,T n n,i,li 
n— 1 i— 1 



(4) 



where n n ^ )CT = <? n i G c n ^ t(J and U > is the repulsive cou- 
pling constant. The electron operators c ni j lCT and c' n i a 
satisfy the periodic boundary conditions for both the leg 
and the rung directions, 



and anticommutation relations, 

{Cm,i,(7i Cn,j,r} = 0, { C m,i,a' c n J.rl = 0- 

The hopping parameters are given by t > 0, Si t 2 = S2,3 = 
s > 0, and S3 f i = /3 s > 0. The strong coupling expansion 
shows that this model at the half-filling case is reduced 
to the Heisenberg model with J\ = At 2 /U , J r = 4s 2 /t/, 
and a = P 2 . 



By performing the suitable unitary transformation, the 
hopping Hamiltonian can be mapped to the following di- 
agonal form: 



^°p = EE E s i(*)4,<,<A*,° 

k i=l 0=1,1 



(5) 



where the wave number k is summed over < k < 2ir. 



The operators dk,i, a an d d\ . a are defined by 

L 3 

7f EE e_,fcno tf^>' ( 6a ) 



n=l j=l 
L 3 



d\,i,<j — ^nj- E E 



n n,j,a > 



(6b) 



n=l j=l 

which satisfy the standard anticommutation relations 

{dk,i,<T,dj j T } = Sk,iSi,jS a ^ T , 

Ka.--'//,,-! {d{ ha ,dl T } = 0. (7) 

The explicit form of the orthogonal matrix is 



O 



where the matrix elements are given by 




1 



u± = — , v± = — , n± = \/2- 
n± n± 



3/3 ± + 8 



(3 2 + 2 ± P^fF+S 
The energy eigenvalues of the one-electron states are 



(8) 



Ei(k) = -/3s + 2tcosk, 
1 



(9a) 



E 2 (k) = -(j3s- s + 8 + At cos ft), (9b) 



E 3 (k) = -{I3s + sy/f3 2 + 8 + 4icos/c). (9c) 

Note that a degeneracy E\{k) = E%(k) appears at /3 = 1 
due to the translational (Sjj — > Sj+ij) and the parity 
(Si j <-> ^j) symmetries along the rung direction. 

For the half-filled case, 3L one-electron states should 
be occupied by electrons with up and down spins. As 
a result, the ground state of the hopping Hamiltonian 
has one, two or three pairs of the Fermi points (kj,kj) 
just on the Fermi sea, depending on the parameters s/t 
and /3. Since the low-energy excitations are given by the 
particle-hole creations around these Fermi points, they 
may be represented by using the Dirac fermions, the left 
mover ipj >a (x) and the right one ipj >a -(x), which are de- 
fined from the electrons around the j-th pair of Fermi 
points. If the j-th. band has no Fermi points in the half- 
filled case, we should neglect ipj^ a (x) and ipj tCr (x). On 
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this understanding, we approximate the original electron 
operators in terms of the Dirac fermions as follows: 



3=1 



{e ik ' x ' a ^, a {x) 



( 10 ) 

where a is the lattice spacing with dimension of length 
and x = an is the continuous position coordinate. 

For this free Dirac fermion system, we take into ac- 
count the effects of the on-site Coulomb interaction ((4]) as 
the perturbation, and we use the non-Abelian bosoniza- 
tion techniques. Following naively the field-theory argu- 
ment in Ref. l30l . one can expect that when the number 
of Fermi-point pairs is odd (even), the spin excitations 
are gapless (gapped) in the half-filled Hubbard tube. In 
particular, in the cases of one or two Fermi-point pairs, 
we can explicitly determine whether or not a spin gap 
exists as follows. 

First, we consider the case of one pair of Fermi points 
k% = and k\ = J. In this case, the interaction (H|) is 
approximated as the sum of an Umklapp interaction and 
two marginal ones 

■Hint ~ J dx gi&i{x) ~ 5262(2:) - .9363(20 + • • 

where 31,2,3 are positive coupling constants proportional 
to U. The Umklapp term is expressed as 

9i(i) =^(i) t ^iW t ^rW^W. ( n ) 

and the marginal interaction between the U(l) charge 
currents is given by 

02 (x) = ipiix)* 1 ipxix^iixY ipi(x), 

where ipi — (^l.t? ^i.i)* ^ is known that the bosonizcd 
form of 01,2 contains only the charge degrees of freedom 
and they open a charge gap when 32 is positive. Then, 
the remaining spin degrees of freedom are described by 
the gapless level- 1 SU(2) WZW theory,±&±&2° This phe- 
nomenon, i.e., the charge-spin separation is well-known 
in the single Hubbard chain model. For this WZW the- 
ory, the third interaction 



©3(2;) = -0i(2;) t |'0i(2;) • ipi{x)^ipi(x), 



(12) 



is known to be marginally irrelevant if 33 > 0. The 
coupling constant 33 is hence renormalized to be zero in 
the low-energy limit. Except for the above interactions 
01 2,3(2;); there is no relevant operator with the invari- 
ance under the one-site translation along the leg, 

1>iA x ) -» e*Vi, CT (z), i>iA*) - e^i, CT (x), (13) 

as in the case of the Heisenberg chain. The spin excita- 
tions, therefore, remain gapless. 

On the other hand, when there exist two pairs of the 
Fermi points, (&i,fci) and (fc2,fc2), the fate of the spin 



excitations is different from the above scenario. In this 
case, the spin sector in the hopping part of the Hubbard 
tube is described by a level- 2 SU(2) WZW theory, which 
is derived from two decoupled Dirac fermionsi^I The 
Coulomb interaction yields several perturbations for this 
theory. For example, applying the non-Abelian bosoniza- 
tion rule^ we find that an interaction, derived from 
Eq. Q, 



^,TW t ^,T(a;)^i,i(i) t ^,i(i) 



(14) 



contains a relevant perturbation for the level-2 WZW 
model, and it is invariant under the translational opera- 
tion 

ipjA x ) -> e zk ^ J: a(x), $j, a (x) -> e %k ^ j>ff (x). (15) 

Particularly for fj = 1, the number of possible relevant 
operators with the translational invariance is increased, 
because of the coincident Fermi points k% — &2 and 
k\ = &2- Therefore, we conclude that any gapless spin 
excitation generally has no chance to survive except par- 
ticularly rare cases [e.g., when the coupling constant of 
Eq. ([14]) is zero]. 
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FIG. 2: Phase diagram obtained from the effective Hubbard 
model ((2]). In the strong-coupling limit (U/t,U/s 3> 1), the 
horizontal axis (3 2 and vertical one {s/t) 2 can be regarded as 
a and J r /Ji, respectively. The effective theory claims that 
the phases (II) and (IV) are gapless and the phase (III) is 
gapful. The phase (I) is possibly gapless. 



Since for the case of three Fermi-point pairs, the in- 
teractions among three Dirac fermions, generated from 
Eq. (j4]), are fairly complicated, it is difficult to analyze 
them and judge whether or not the spin excitation can 
survive as gapless. However, as we already stated, it can 
be expected from Ref. [3(3 that the spin excitation is pre- 
sumably gapless in this case of three Fermi-point pairs. 

From these arguments, we can draw the ground-state 
phase diagram of the half-filled Hubbard tube as shown in 
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Fig. [2j The phase (I) has three pairs of the Fermi points, 
the phases (II) and (IV) have one pair, and the central 
phase (III) has two pairs. Therefore, we can predict that 
the phases (II) and (IV) have gapless spin excitations, 
whereas the phase (III) possesses a spin gap. If {s/t) 2 
and 1 are, respectively, replaced with the coupling ratio 
Ji/J r and the asymmetric parameter a in Fig. [51 one 
may interpret the phase diagram as that of the S = h 
three-leg spin tube (JTJ) . Although we have treated the on- 
site interaction Q perturbatively, we have assumed that 
the weakly-interacting case (U/t,U/s <C 1) is smoothly 
linked to the strongly-interacting one, which is nothing 
but the spin tube (P). 

We here note that the gapful phase (III) is predicted 
to be extended around the line [3=1 for a finite s/t. 
In the limit s/t — > oo, the both the left- and right-side 
phase boundaries of the region (III) converge to (3 = 1 . 
This narrowing of the phase (III) is consistent with the 
numerical result o 11 ' 12 in the strong-rung-coupling limit 
J r /Ji — ► oo. We will discuss this limit in more detail in 

Sec. Unci 

Finally, we briefly argue the universality classes of 
the phase transitions at two phase boundaries, (II)- (III) 
and (III)- (IV). For the level- 1 SU(2) WZW model in 
the phases (II) and (IV), the most relevant perturba- 
tion is the marginal current-current interaction ; 18 ' 19 ' 20 
Eq. (|12p . which is invariant under the translational and 
spin-rotational operations. Since it possibly becomes 
marginally "relevant" when parameters are finely tuned 
and then g$ becomes negative, we speculate that the 
transition from the phases (II) or (IV) to (III) is caused 
by this marginal term. Therefore, the transitions are 
expected to be in the BKT universality class. This spec- 
ulation may be naturally accepted if we recall the follow- 
ing two known results of the S = | zigzag Heisenberg 
chain, namely, the spin chain with the nearest- and the 
next-nearest-neighbor interactions; 32 (i) When the next- 
nearest-neighbor hopping is sufficiently small in the half- 
filled electron system on zigzag lattice, one obtains one 
pair of Fermi points. The spin excitations are therefore 
described by a level- 1 SU(2) WZW model, like the phases 
(II) and (IV). (ii) It has been numerically shown 32 that 
when the next-nearest-neighbor interaction is increased 
in the zigzag spin chain, the BKT transition takes place 
and the ground state changes into a dimerized state from 
a Tomonaga-Luttinger liquid described by a level-1 SU(2) 
WZW model. 

On the basis of this speculation, we numerically an- 
alyze the phase transitions between (II)- (III) and (III)- 
(IV) in Sec. [TV] 



B. Gapless phase for J r <C J[ <C Ji 



In Sec. IIII A[ we have obtained a qualitative phase di- 
agram of the spin tube ((T]) as shown in Fig. [51 However, 
some subtle points still remains, and particularly it is 
doubtful whether or not there is the gapless phase (I). 



To partially resolve these issues, we here focus on the 
extreme situation J r <C J[ < Ji, (a = J' r /J t ~> 1), 
which corresponds to the right-lower regime in Fig. [5] In 
this regime, we prove the existence of the gapless phase 
(I) and discuss the phase transition between the gapless 
phase (I) and the gapful phase (III). 

Before embarking on our analysis, we sketch the sce- 
nario of this subsection. We introduce three level-1 
SU(2) WZW theories for the three decoupled Heisenberg 
chains ; 18 ' 19 ' 20 and then we treat their rung couplings as 
the perturbation, because the weak-rung-coupling regime 
J r , J[ -C J\ is considered now. The first and the third 
chains are coupled to each other with J r ' which is much 
stronger than two remaining couplings J r . It is well- 
known that the J[ coupling involves a relevant interac- 
tion with conformal dimensions (4, i) in the two coupled 
WZW models [see Eq. ([20]) ], and it produces an energy 
gap^ 3 - On the other hand, like the case of one Fermi-point 
pair in Sec. IIII A[ the WZW model for the second chain 
has the marginal irrelevant current-current interaction, 

A 2 $ 2 , (16) 

with a finite negative coupling constant A2 < 0. The op- 
erator $2 is equivalent to Eq. (fl5|) . if we use the Dirac 
fermions (V'i.ct, V'i.o-) to describe the second chain. The 
negative sign makes Eq. (|16[) irrelevant and the second 
chain is gapless. A weak rung coupling J r between this 
WZW model and the massive theory for the two coupled 
chains must give a correction to the coupling constant A2 . 
If J r is sufficiently small, the sign of A2 < would not 
change and the gapless excitation is preserved. These ar- 
guments convince us that the gapless phase (I) definitely 
exists. Furthermore, we expect the phase transition from 
the gapless phase (I) to the gapful phase (III). Namely, 
if J r exceeds a critical value, the coupling constant A2 
might change to be positive. In this case, the marginal 
operator (|16p becomes relevant, which produces an exci- 
tation gap. This transition between the phases (I) and 
(III) is nothing but of the BKT type ; 24 ' 25 as in the zigzag 
Heisenberg chain^ 

Now, let us calculate the correction to A2 in an ex- 
plicit way to confirm the above sketch. We represent 
the partition function of the three coupled WZW mod- 
els through the path integral formalism. Note that 
the SU(2) spin-rotational symmetry is conserved auto- 
matically in the language of the WZW models, i.e., 
the non-Abelian bosonization. The three WZW models 
can be represented in six Dirac fermions ^>i_ a (z,z) [i — 
1,2,3, and a =T>4] and three ghost bosons 4>i(z,z) [i = 
1, 2, 3]. Here, z = vt + ix and z = vt — ix are the chiral 
coordinates (t: imaginary time) and the index i denotes 
the chain number of the spin tube {!]). We employ this 
free-field represent at ion j 34 1 3 5 1 36 The Lagrangian density 
for the decoupled three chains is given by 



£ Q (z,z) = 2j2 



E 



(17)" 



6 



where we have decomposed the Dirac fermions into the 
linear combination of the left mover ^ ia (z) and the 
right mover *i, a (z) as ^ i>a {z,z) = Vi, a (z) + * 4 , a (z). 37 
The nonvanishing two-point correlation functions of these 
fields, calculated from this free Lagrangian (fP7|) . are 



( e **«{*.*) e -»* ) = S tJ \z-w\. 



z — w 



(18a) 

(18b) 
(18c) 



The anomalous correlation of the ghost bosons should 
be noted. These ghost fields kill the charge degrees of 
freedom and extract the gapless spin degrees of freedom 
in the Dirac fermions. The primary field g l a p{z,z) with 
the conformal dimension (4, i) in the i-th level-1 SU(2) 
WZW theory is represented in terms of these free fields: 

gi {z^=% >a {z^%^y MZ ^- 

Furthermore, the spin operators are represented as 

5i,„/aw Ji{z,z) + (-l) n Ni(z,z), 

in which the smooth and the staggered parts are given 
by 



Ji(z,z) = tftf^ + fctf^ 



(19a) 



Ni(z,z) = C (tft-^e** +h.c.J , (19b) 

with ^S>i = '(^i.-f, Here Co is a nonuniversal con- 

stant. Applying Eqs. ()18a|) - (|18c|) . one can evaluate the 
asymptotic forms of several correlation functions for the 
field g l a p and the spins Si >n . Among the spin-rotational- 
and the translational-symmetric operators, the most rele- 
vant coupling between the i-th and the j-th WZW models 
is 



$ij(z,z) = Ni(z,z) -Njfaz), 



(20) 



with the conformal dimensions (o> 5)- As we already 
mentioned, this relevant term $31 (z, z) produces the en- 
ergy gap in the two strongly coupled WZW models. 
The WZW model for the second chain interacts with 
these coupled WZW models via the weak rung coupling 
^12(2, z) and $23(2, z). The explicit form of the marginal 
operator (TlB|) of the second WZW model is 



(21) 



From these materials, the total Lagrangian density for 
the asymmetric spin tube under the condition J r <C J' r <C 
Ji is written as 

C(z,z) = C Q (z, z) + A 3 i$ 3 i(z,z) + \i2$ 12 (z,z) 

+X 2 3^23(z,z) + X 2 <i>2(z,z) + --- , (22) 



where the coupling constants are A31 oc J' T , A12 = A23 c>c 
J r , and A2 < 0. Therefore, the partition function is writ- 
ten in the path integral over Grassmann and boson fields 

as 



x exp 



fz 
~2n 



C(z,z) 



Under the condition J r <C J[ <C Ji, the low-energy 
physics must be governed by the second chain weakly 
coupled to the other two chains. To obtain its effective 
theory, we may integrate out the massive degrees of free- 
dom ^i ta (z,z) and cf>i(z,z) with i = 1,3. To calculate 
the correction to the coupling constant A2, we regard the 
relevant term $31(2,2) as an unperturbed Lagrangian 
and expand the partition function in all the other oper- 
ators Co, $12(2, z), $23 (z, z), • • • . This expansion can be 
performed by a lattice regularization of two-dimensional 
Euclidean spacetime. In this regularization, the partition 
function is represented as the following multiple integra- 
tion over countable variables: 



Z 



n n n {a a d*u*»)<**u*o 



(23) 



d* i , Q (z„)d* l t Q (z„)# i ( 



Z7T 



where I = {1,2,3}, S = {T,|}, and N = {(m,n 2 ) £ 
Z 2 1 1 < ni, ri2 < L} is a finite set of integer pairs 
with a large number of elements. For an integer pair 
n = (^1,712) € N, we define discretized coordinates 
(z„,z„) = a(ni + iri2,ni — in 2 ) for (z,z) with a lattice 
spacing parameter a. For a finite number of Grassmann 
variables, the Taylor expansion is reduced to a finite sum- 
mation, 



exp 



£ \%n ; Znj 
Z7T 



K 



1 



^ k\ 

k=l 



£(z n: z n ) 

Z7T 



with a certain positive integer K < 12 jiV j , due to the 
nilpotcncy of the Grassmann variables, 

^iA Z nf = 0, ^i, a (Zn) 2 = 0, 

*Uv> 2 = °> *U^) 2 = o. 

The following integration formula of Grassmann vari- 
ables, 



n n n {^^^A^d^A^^^A^d^iA^. 
n n n {^lA^iA^A^iA^)) 

n£M iG J o£T 

1 : M — N, J = I,T = S 
: otherwise, 



it 



(24) 
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is important to calculate this expansion. In addition, 
using an equality for two Grassmann variables 

and a trace formula of the Pauli matrices 

i._r_<l- r 2 t J„2l \ " a 2 b J2 0X ab 

one can reduce products of Pauli matrices and Grass- 
mann variables to the following single term, 

= n*lT*U**.T*U*l T *U**.T*U?-V' (25) 



These formulas (1241) and (|25[) make it easier to calculate 
the expansion of the partition function (|2"3")l . As a result, 
we obtain 



rt67VaeS 
d*2,a(-2ra)^2 a (Zn)d4> 2 (z n , Z n )] 



n{- 



2tt 



x I I i exp 

r 15 A 31 5 A3JA12A23 



$2(z«,Zn)H } 



L32(2tt) 4 16 (2tt) 5 



exp 



2n \ ' 3A 31 



• (26) 



Here, we have neglected several terms of charge degrees 
of freedom. The final expression in Eq. (|2"6")l clearly in- 
dicates that the correction to A2 is 2Ai2A23/(3A3i) = 
2J r /(3a), and the phase boundary between (I) and (III) 
is given by 



J r = --A 2 a. (A 2 < 0). 



(27) 



Thus, the gapless excitation in the second chain is pre- 
served and the gapless phase is expanded under the con- 
dition J r <C J' t <C J\. Namely, we have proved that 
the phase (I) is exactly present at least in the region 
J r <C Jr <C Ji- Furthermore, we find that if J r is large 
enough to change the sign of A 2 + $2(2, z) becomes 
relevant and an energy gap appears. Since this phase 
transition is induced by the marginal operator <& 2 (z, z), 
it belongs to the BKT universality class, as in the S = | 
zigzag Heisenberg chain. This gapped state must corre- 
spond to the phase (III) in Fig. [21 In the gapful region 
under J r <C J r ' <C Ji, the second chain must be dimer- 
ized. On the other hand, as we have often mentioned, 



the gapful phase in the vicinity of the symmetric line 
a = 1 has the valence-bond order. Since both the orders 
break the same translational symmetry, the state with 
the dimerized second chain would smoothly change into 
the valence-bond ordered state when we vary a from a 
large value to unity in the gapful phase (III). 

Finally, we notice that in the prediction l|27p . the gap- 
less phase (I) becomes wider for larger a, contrary to 
the phase diagram, Fig. [2] Two reasons for this contra- 
diction are immediately found. First, we should not pre- 
cisely trust the location of the phase boundaries in Fig. [21 
which are drawn just by counting the number of Fermi 
points. The argument in Sec. IIII A1 is only qualitatively 
correct to determine the phase boundaries. Particularly, 
since the gapful phase (III) sandwiched by two gapless 
phases (I) and (IV) is quite narrow for large a in Fig. [21 
then we cannot claim the existence of the gapful phase 
(III) for a ^> 1 by using Fig. [21 Therefore, the phases 
(I) and (IV) might be smoothly connected in the large-a 
region. For such a case, the transition between two re- 
gions (III) and (IV) is expected to be also of the BKT 
type. Second, the three Dirac fermions in the region (I) 
do not always imply the gapless phase ) 27 i 28 as discussed 
in Sec. IIII Al The coupled three Dirac fermions possibly 
produce an energy gap by the strong frustration, if a ap- 
proaches the symmetric point a = 1 . The result ([2"T]) also 
shows this tendency. 



C. Energy gap for J r 3> Ji and J r ~ J' r 

In this subsection, we discuss the existence of an energy 
gap for a sufficiently strong rung coupling J r ^> J\ and a 
sufficiently weak asymmetry \a— 1| <C 1 except at J\ = 0. 
Namely, we explain that the phase (III) has a finite width 
along the line J r / Ji = const 3> 1. 

As we mentioned before, Kawano and Takahashi 6 have 
obtained the spin-orbital type model as an effective the- 
ory for the symmetric (a = 1) spin tube with a strong 
rung coupling. The Hamiltonian is written as 



rrsyro 
n cff 



3J r 



3+1 



l+4(r. 



J T 3 + l~ 



~ T 3 T l + l' 



. (28) 



3 = 1 



which is the result of the first-order expansion in Ji/J r . 
The orbital Pauli matrices t" (y = x, y, z) are defined for 
the states with respect to the left- and right- handed spin 
configurations on each rung (see Refs. la and IT3). From 
this theory (|28|) , Kawano and Takahashi have shown that 
the ground state has a valence-bond order with the trans- 
lational symmetry (Si : j — > Si j+i) spontaneously broken, 
and a finite excitation gap exists. It is believed that the 
energy gap is generated by the strong coupling between 
the spin and orbital degrees of freedom in the effective 
Hamiltonian ([28]) . Obviously, either of them becomes 
gapless, if the other is frozen. 

Quite recently, Nishimoto and Arikawa^ have ex- 
tended the effective theory to the asymmetric case as 
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follows: 

ff e T m =^:r+(«-i)E^- (29) 

For simplicity, we have re-defined the orbital Pauli ma- 
trices t" through a unitary transformation. Remarkably, 
the asymmetry induces a transverse external field cou- 
pled to the orbital spins. A sufficiently strong asymme- 
try | a — 1| 3> 1 hence yields the saturated polarization 
of the orbital spins. In this case, the effective Hamil- 
tonian (|29| is reduced to that of the antiferromagnetic 
Heisenberg chain. The spin degree thus survives as a gap- 
less excitation. Next, let us focus on the region around 
the symmetric line a ~ 1. In the limit J\/J r — ► 0, the 
ground state has a saturated orbital spin at any asym- 
metric point a 7^ 1. For a finite J±/J T , however, we 
expect an extended gapful phase around a = 1, if the en- 
ergy gap exists at a = 1 due to the coupling between the 
spin and the orbital degrees of freedom. The finite en- 
ergy gap does not vanish by an infinitesimal external field 
a — 1. In other words, the magnetization process of the 
orbital spins should show a zero-magnetization plateau. 
Therefore an energy gap would also be present for suffi- 
ciently weak asymmetry \a — 1| <C 1. 

From these arguments, we find that the gapful phase 
surely exists in a finite region under the conditions 
J r ^> J i and J? ~ J r ', and its width becomes smaller 
with decreasing J\/J T . This is consistent with the result 
obtained in Sec. IIII Al 



IV. NUMERICAL ANALYSIS 

In this section, we numerically analyze the asymmet- 
ric spin tube |T]), taking into account the results in the 
preceding Sec. 11111 We show that the arguments based 
on the effective theories are in good agreement with the 
results of the numerical calculations. We will use results 
of the numerical diagonalization up to L — 10 for the 
periodic boundary system, and those of the DMRG up 
to L — 128 for the open boundary system. 

A. Phenomenological renormalization 

A useful order parameter to determine the phase 
boundaries between the gapless and the gapful phases in 
Fig. [21 is the spin gap A, which is the energy gap between 
the singlet ground state and the triplet excited state in 
the finite, but large system. We calculate it by means of 
the DMRG up to L = 128, where we do not see such a 
significant open boundary effect as local edge excitation. 
In actual DMRG computation, the number of retained 
bases is up to m — 300, within which well convergence 
is achieved for the scaled gap. In Fig. [3] (a), the scaled 
gap LA is plotted versus a for Ji = 0.5. It indicates 



that the spin gap is just open in a tiny region a ~ 1 and 
it rapidly vanishes away from a = 1. Thus we find two 
critical points a c \ and a C 2 (a c i < 1 < a. C 2), which are 
expected to be of the BKT type because of the wide gap- 
less regions outside the gapful phase, and the discussion 
in Sec. QUI 

In order to determine a phase boundary of the usual 
second order phase transition, the phenomenological 
renormalization equation LiA^^ac) = L2&-i J2 (a c ) is of- 
ten used effectively. For the present critical point a c , 
however, this type of phenomenological renormalization 
has no clear crossing point. It seems from Fig. [3] (a) that 
the scaled gap LA increases with increasing L in both 
gapless and gapful phases. This is because the scaled 
gap is an increasing function with respect to L not only 
in the gapful phase but also in the gapless phase, since 
the finite-size gap must have the logarithmic size cor- 
rection term ~ — 1/logL. Here, we should recall that 
the logarithmic correction normally vanishes just at a c 
due to the SU(2) symmetry in the c = 1 CFTM There- 
fore, instead of using the crossing point of the scaled 
gaps, we can estimate a c as a point where the size cor- 
rection is minimized. The difference of the scaled gap 
between two system sizes L\ and L2 is plotted versus a 
for Ji = 0.5 in Fig.[3](b). The minimum of the difference 
L\Al 1 — L2Al 2 has a very small L dependence. Note 
that the minimum value L±Ai J1 —L^A^ decreases as the 
size increases. This phenomenon and the assumption of 
the BKT type transition suggest that this minimal value 
approaches zero as the system size increases. This is quite 
reasonable if we suppose the most important finite-size 
correction to the scaled gap LA next to 1/ log L term is 
order of 1/£ 2 j 32 i 38 i 39 We thus determine a c from the min- 
ima for two large systems with L\ — 96 and L2 = 128 for 
J 1 < 2. The estimated a c i and a C 2 are shown as crosses 
in Fig. 2] They correspond to the phase boundary be- 
tween two regions (II) and (III) and that between (III) 
and (IV), respectively. At least these boundaries for the 
strong-rung-coupling regime J\ < 2 are precise enough 
to justify that a finite gapful phase (III) exists. However, 
it is difficult to obtain a c for J\ > 2, because the DMRG 
calculation is not well converged there. 

In order to determine the phase boundaries for the 
weak-rung-coupling regime J\ > 2, we use the mini- 
mum points of L\Al 1 — L2Al 2 calculated by the nu- 
merical diagonalization up to L = 10 under the periodic 
boundary condition. The estimated phase boundaries for 
(L1.L2) — (6,8) and (8,10) are plotted as long-dashed 
and dashed curves, respectively, in Fig. 3J In addition, 
the infinite-!/ curves extrapolated assuming the size cor- 
rection is proportional to 1/L in both directions of J\ 
and a are also shown as solid curves in Fig. [5] At least the 
phase boundaries (II)- (III) and (III)- (IV) are consistent 
with the DMRG results for J x < 1. The boundary (III)- 
(IV) is, however, significantly deviated from the DMRG 
estimation for J\ *~ 1. This discrepancy is supposed to 
be due to the error of extrapolation. This analysis also 
justifies the existence of the phase (I). However, the error 
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FIG. 3: (a) Scaled gap calculated by DMRG for Ji = 0.5. (b) 
Difference of the scaled gaps between two systems with sizes 
L\ and L2. 
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FIG. 4: Ground-state phase diagram of the asymmetric three- 
leg spin tube (TTJ, derived from the numerical analysis. The 
phases (I)- (IV) correspond to those in Fig. [2] The cross points 
are determined by the DMRG. The three kinds of lines are 
done by the numerical diagonalization. The circle points are 
obtained from the level-spectroscopy method in Sec. IIV CI 



nents 77 and er t appear in the spin correlation functio: 



,36,40 



of extrapolation becomes larger as we approach the line 
1/Ji = in the case of a < 1. Thus it is difficult to 
conclude that the phase (I) really exists for a < 1 within 
the present numerical demonstration. The boundary (I)- 
(III) will be discussed later. It is also difficult to confirm 
the boundary (I)- (II), and the phase (I) might combine 
with the phase (II) in a certain regime with a < 1. 



B. Conformal field theory analysis 

The numerical analysis based on CFT is generally effi- 
cient in investigating the feature of quantum phase tran- 
sitions in 1+1 dimensional systems. The conformal in- 
variance gives the L dependences of the ground-state en- 
ergy .Bo and the energy gaps for triplet excitations A t 
and that for singlet excitations A s as the forms, 

Eq ttv s c 

T °° _ IXH ' (30a) 
*-?(--SiE + -)- ,3 ° b) 

?("-&+■■■)• ^ 

where v s is the spin wave velocity, c is the central charge, 
and 77, <7t and c s are the critical exponents. These expo- 



(StfiStj) - (-l) J (log|j|r|jr, {v = x,y,z.), 

where \j\ ^ 1. As discussed in Sec. lIIIl the gapless phases 
(II) and (IV) are predicted to be described by the level- 1 
SU(2) WZW model. This model indicates the universal 
constants c = 1 , rj = 1, <7 t = i and ct s = — | . 

The exponent 77 can be determined from the triplet 
and singlet excitation gaps in the finite size system, by 
using the relations rj — LA t /Trv s +- ■ ■ and 77 = LA s /7tt; s -|- 
• • • . Using the results from the numerical diagonalization 
for L = 8 and 10, we estimate c, LA t /Trv s = rj t and 
LA s /nv s = 7] s independently, shown in Fig. [5] for J\ = 
0.3. Here, we have evaluated the spin wave velocity by 
v s = (Ex — E )/ki, where Ex is the lowest energy of the 
eigenstate with the smallest nonzero wave number k\ — 
=r-. The velocity is generally nonuniversal and depends 
on the couplings J% and J T . 

In Fig.0 we also depict a special average (3L At /ttv s + 
LA s /ttv s )/A of the singlet and the triplet gaps such that 
the dominant finite-size logarithmic corrections cancel 
out each other.— This average and c seem to be al- 
most unity in the gapless phases in Fig. [5] as observed 
in the S = ^ antiferromagnetic Heisenberg chain with 
the next-nearest-neighbor interaction.— These results are 
completely consistent with the expected BKT transition. 
We therefore conclude that both gapless phases (II) and 
(IV) are governed by the WZW model, and the transition 
between (II)- (III) and that between (III)-(IV) belong to 
the BKT universality class. 
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FIG. 5: Central charge c and critical exponents LAt/nv s and 
LAs/irvs including finite size corrections for Ji = 0.3 . These 
are estimated from the numerical diagonalization up to L — 
10, based on the size dependence of low-lying energy spectra 
predicted by the CFT. 



Now, we note the level crossing between L A t / nv s and 
LA s /irv s around a = 1 in Fig. [5] This level crossing 
implies the appearance of another singlet ground state 
as a reflection of the valence-bond order in the thermo- 
dynamic limit. This is a clear numerical evidence of the 
extended gapful phase (III) predicted by several effective 
theories. On the contrary, Nishimoto and Arikawa 12 have 
claimed that the gapless phase is extended everywhere 
except at the point a = 1, by means of the DMRG anal- 
ysis. Our observation does not agree with their claim. 
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FIG. 6: Critical points a c evaluated by the level-spectroscopy 
analysis for Ji = 0.2, assuming the size correction is propor- 
tional to 1/L 2 . 




FIG. 7: Critical value of Ji for the regular triangle spin tube 
(a = 1) evaluated by the level-spectroscopy analysis. 



C. Level-spectroscopy method 

The level-spectroscopy metho d 32 i 42 i 43 ' 44 is a very pow- 
erful tool to determine the critical point of the BKT 
transition in one-dimensional quantum systems. For the 
SU(2)-symmetric cases including the present spin tube, 
its strategy becomes easier £1 According to this method, 
the critical point can be determined as an intersection be- 
tween the singlet and the triplet excitation gaps, where 
their logarithmic finite-size corrections vanish. The phase 
boundaries (II)- (III) and (III)- (IV) estimated by this 
method are shown in Fig. [BJ Here we have applied the 
results of the numerical diagonalization up to L = 10 for 
Ji = 0.2. The numerical data are well converged to those 
of the thermodynamic limit by use of the 1/L 2 extrapo- 
lation. The 1/L 2 extrapolation is justified by considering 
the most important finite-size correction to the excitation 
gaps A next to the logarithmic termi^ 2 - 

Several critical points estimated by the level spec- 
troscopy are plotted by open circles in Fig. 2J We find 
that at least for Ji <C 1, the results are in good agree- 
ment with the phase boundaries (II)- (III) and (III)- (IV) 
evaluated in the preceding Sec. I IV Al 

On the other hand, when we consider the transition 
between two phases (III) and (I) along the line a — 1, 



the finite-size correction to the critical point is too large 
to determine the precise value of the infinite-length limit, 
as shown in Fig. [7j The extrapolated value of the critical 
Ji results in 1/Ji c = 0.51±0.45. It is therefore difficult to 
conclude whether J lc is finite or zero for a = 1, namely, 
whether or not the gapless phase (I) still survives for 
small a. [we have already shown the presence of the 
phase (I) for large a 3> 1 in Sec. IIII 51 ] 



V. CONCLUSIONS AND DISCUSSIONS 

We have studied the ground-state phase diagram and 
the quantum phase transitions of the S = i three-leg 
asymmetric spin tube models, defined by Eq. JI]). In 
Sec. IIII Al based on the Hubbard model on the tube 
lattice, we have proposed an effective theory to draw 
a global phase diagram, Fig. [21 in the parameter space 
(Ji,a). Three gapless phases (I), (II), and (IV) and one 
gapful phase (III) are found by counting the Fermi points. 
This effective theory indicates that the level-1 SU(2) 
WZW model describes two extended gapless phases (II) 
and (IV) which are separated by a extended gapful phase 
(III) around a = 1. In Sec. IIII B[ applying another an- 
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alytical strategy based on the non-Abelian bosonization, 
we have proved that the predicted gapless phase (I) is 
exactly present at least for weak rung couplings J\ ^> J r 
with a strong asymmetry a»l. Furthermore, we have 
argued that the gapful region (III) is surely extended 
(although narrow) for the case J r 3> Ji and J r ~ J r ' in 

sec. mrn 

Following these results of effective theories in Sec. IIII1 
we have numerically analyzed the quantum phase transi- 
tions of the spin tube ([TJ) in Sec. IIV1 The phenomenolog- 
ical renormalization approach based on the DMRG and 
the numerical diagonalization has enabled us to draw the 
global phase diagram. The numerical results are quali- 
tatively in agreement with those of effective theories. In 
addition, we have raised the validity of the phase diagram 
by means of the numerical finite-size scaling arguments 
based on the c = 1 CFT. We have confirmed that the 
phase transitions (II)- (III) and (III)- (IV) belong to the 
BKT universality class quantitatively. Here, we make a 
comment on that the gapped region estimated in Ref. Il2l . 
where an unconventional power law fitting with respect 
to L is employed, is well inside of the present phase III. 
We should, however, recall that the gap near the BKT 
transition is exponentially small and thus the gap is dif- 
ficult to be detected by the phenomenological renormal- 
ization approach, which usually overestimates the gapless 
region. This suggests that the gapped region in Rcf. 12 
becomes smaller than the present phase diagram (Fig. [4]), 
although the DMRG data itself may be consistent with 
each other. 

The semi-quantitative phase diagram of the spin 



tube ([T]) has been constructed in this study (see Figs. [2] 
and 2]), but some subtle issues are still remaining, e.g., 
(A) the topology of the phase boundaries along the 
two lines a = and 1/a = 0, (B) how widely the 
gapless phase is extended, etc. In order to resolve 
these problems, more sophisticated approaches would be 
necessary. 
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